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Abstract. We present some elliptic beta integrals with a base parameter on 
the unit circle, together with their basic degenerations. 



1. Introduction 

The theory of generahzed gamma functions has been set up by Barnes |Ba| . 
A slightly different approach was advocated by Jackson, who considered the basic 
gamma function depending on one base parameter q and the elliptic gamma function 
depending (symmetrically) on two bases p and q p]. For a long time, only the 
first of these generalized gamma functions was appreciated in the literature |AAR,j . 
Recently, however, the elliptic gamma function also got appropriate attention after 
the work of Ruijsenaars |Rul| . who introduced it in the context of intcgrable systems 
and investigated some of its properties. A further study of the function in question 
was conducted by Felder and Varchenko [FVj . A modified elliptic gamma function, 
which admits analytic continuation in one of the base parameters, e.g. q, onto the 
unit circle |g| = 1, has been introduced recently by one of us in S2 . 

In this paper we study beta type integrals on the unit circle \q\ = 1 built of 
modified elliptic gamma functions, as well as their basic degenerations. The first 
exact beta type integration formula involving the conventional elliptic gamma func- 
tion was discovered in SI . Various multidimensional generalizations of this elliptic 
beta integral associated with the Cn and An root systems have been investigated 
in |DS1I inS2l rni IS2 | . A general theory of theta hypergeometric integrals on tori 
and the Jacobi theta function generalizations of the Meijer function was developed 
in |S2| . The beta integrals considered below should be thought of as the |g| = 1 
counterparts of the integrals in |Slj and |DS1| . Recently, the Bailey's technique of 
deriving identities for series of hypergeometric type |AAR,| has been generalized to 
the level of integrals jS8j . It can also be extended to the |g| = 1 integrals under 
discussion. 

2. Preliminaries: 

THE Jacobi theta function and the elliptic gamma function 

The main underlying structural object of this paper is a Jacobi type theta func- 
tion defined as 

oo 

9{z;p) = {z;p)oo{pz~^;p)oc, {a;p)oo = Y[{1 - ap"), (1) 

with z,p £ C, \p\ < 1. It satisfies the transformation properties 

eipz;p)^eiz-';p) = -z-'e{z;p). (2) 
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Evidently, 0{z;p) = for z = p"\ m e Z, and 0{z;Q) ~ 1 — z. If we denote 
p = e^^^'^ , Im(T) > 0, then the standard Jacobi 6'i-function |WW| is related to 

9{z;p) as 

oo 
n— — oo 

= pi/8,e— n(p.p)^0(e2.».p)^ (3) 
The complex function 9i{u\t) is entire, odd, and doubly quasiperiodic in u 

9,iu+l\T) = -9,{u\T), 

9i{u + t|t) = -e-"'^-2^™0i(u|T). (4) 
The modular PSL{2, Z)-group, 

ar + b u 



CT + CT + d' 



(5) 



with a,b, c, d G 7j and ad^bc = 1, is generated by the two transformations r t+1, 
u —>■ u and r — > —t~^, u —>■ ut~^ . Its action on the Jacobi thcta function is 
determined by 



i(u|t + 1) = e"/40i(w|r), (6a) 
i(-| - i) = -i(-*T)i/2e™'/^0i(u|T). (6b) 



(Throughout this paper the sign of the square root is fixed in accordance with the 
principal branch with the cut chosen on the negative real axis.) From the second 
of these relations, combined with the modular transformation law for the Dedekind 
?]- function 

e (^e ' ;e = (^*^) e 12 (e ;e , (7) 

one readily deduces a corresponding modular transformation formula for the 9(z;p) 
function 

9{e^'''^;e-^'''i) = -ie^*('^+^+B^+^-")6i(e2™; e^'^'^). (8) 

The elliptic gamma function T{z;q,p), \q\, \p\ < 1, is related to the above theta 
function through the difference equations 

^{qz;q,p) ^ 9{z;p)r{z;q,p), r{pz; q,p) = 9{z; q)T{z; q,p). (9) 

It is given by the explicit product representation | Kul| 

nz;q,p)= n 1 . r - (10) 
1 — zqip'^ 

j,k=0 

From this representation the following reflection relation is immediate 

^{z;q,p)Tiz-^;q,p) = — -J", _^ (11) 

9{z;p)9{z i;g) 
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3. The modified elliptic gamma function 



The modified elliptic gamma function introduced in jS2| is constructed as a 
product of two elliptic gamma functions of the form in IjlUII . corresponding to two 
different choices of bases. It is convenient to pass to an additive formulation by 
introducing three pairwise incommensurate quasiperiods coi, L02, 0J3, and write 

27ri^ 27ri^ 

q = e ^2 , p — e , r — e "i, 

q — e "1, p — e , r = e "3, (12) 

(i.e., T ~ 0^3/1^2) • The tilded bases q,p, and f are the respective modular transfor- 
mations of and r. For Im(a;i/w2), Im(a;3/wi), Im(a;3/w2) > (so \q\, \p\, \r\ < 
1), the modified elliptic gamma function is now defined as S2 

00 — 27rz-^ 7 + 1 k-\-^\/^ 27Ti— — 74-1 h\ 

Giu;u.) = n P (13) 

(l-e -2qJp'^)(l-e -1 gj^'^+i) 

It satisfies three difference equations 

G{u + uji;uj) = 0{e^'^'^;p)G{u;uj), (14a) 
G(u + W2;'^) = e(e'"'^;r)G(w;a;), (14b) 

Giu + ..3; u;) = ^J\^j7:'^\ Giu; u:), (14c) 

the latter of which can be rewritten with the aid of modular transformation JHl as 
G{u + UJ3- Lj) = e-"-^=^'^("^"')G(u; w), (14d) 

where 

-U^ U U UJi UJ2 I 

B2,2{u;uj) = + + + 

UJ1UJ2 OJi LO2 bUJ2 DOJi 2 

Such a system of three difference equations determines the meromorphic function 
G{u](jj) up to a multiplicative constant (which is the only meromorphic function 
with three pairwise incommensurate periods a;i,2,3). Similar to the r(z;g,p) func- 
tion, the G(u; uj) function satisfies a simple reflection equation given by 

G(w;a;)G(-M;a;) = 5—, ^— . (15) 

0{e "i;r)t^(e '•'2;p) 

If we fix the quasiperiods wi, UJ2 such that Im(a-'i/ci;2) > 0, and take UJ3 to infinity 
in such a way that Im(u;3/a'i), Im(a'3/a'2) — > +00 (so p,r ^ 0), then we obtain 

2 (^2-iTiu/u2. 

lim — = Siu] 0^1, CJ2) — , n ■ I — I ■ (16) 

P,r^oG(u;u;) v ' ^^2mu/u>rq.q^^ ^ > 



In the modern time, the function S{u;u:) was introduced by Shintani |Sh| . It is 
related to the Barnes double gamma function and is called the double sine function 
|Ku| . Its various properties and applications are discussed, e.g., in [F ||JM[|RLS|IM1 
INUIIRulllRu2| . Some g-beta integrals expressed in terms of S{u; uj) were considered 
in |F!r^ lPTll^lT|. 

One of the main properties of the double sine function consists of the fact that 
it can be extended continuously in the quasiperiods uji^2 from the upper half plane 
Im(wi/aj2) > (so |g| < 1) to the positive real axis UJ1/CU2 > (so |g| = 1), 
the resulting function still being meromorphic in u. A similar situation holds for 
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the function G{u;iv), as can be verified by expressing it in terms of the Barnes' 
muhiple gamma function of the third order |S2| . The following theorem provides 
a convenient representation of G{u; ui) detailing explicitly its structure when q lies 
on the unit circle. 

Theorem 1. Let Im[ijj^/uji), Im{uj^/uj2) > and Im[ijJi/L02) > 0. The analytic 
continuation of G{u;uj) Ijl4dp from the domain Im{u!i/uJ2) > to the boundary 
LU1/UJ2 > is given by the meromorphic function 

Giu;u:) = e^"'^(")r(e-'"^; f (17a) 

where P{u) is the following polynomial of the third degree 

««> ^ 3^ (" - 1 1 («^ - " E - + ^ t ^ (I'M 

\ n— 1 / \ n— 1 n— 1 / 

Proof. Let us first assume temporarily that Im(ct;i/a;2) > 0. We denote the right- 
hand side of (|17a|l by f{u). It is easy to see that 

^ e-(P(u)-P(«+-,))0(e-2..^.^^ ^ 0(e2-^;p), (18) 

/(") 

as a consequence of the modular transformation for theta functions in (jS]). The 
function f{u) thus satisfies equation H14a(l . By symmetry, it satisfies (|14b|l as 
well. Analogously, we have that f{u + lo^)/ f{u) — e'^*(^(")~^("+"3)), which is seen 
to coincide with Ijl4d|l . The independence of the quasiperiods wi,2,3 over Q now 
implies that G{u;uj)/ f{u) must be constant. Its value is equal to one, because for 
u — {uji + LU2 + ^3)72 we have that G{u; lj) — f{u) = 1. The theorem then follows 
upon analytic continuation of the right-hand side of H17a|l to the positive real axis 

UJi/u}2 > 0- □ 

Remark 1 . In |FV| , modular transformation properties of the elliptic gamma func- 
tion were investigated. For |q|, |r| < 1 the equality in H17a(l follows from one of 
these modular transformations. Function T{z]q,p) has a pointwise limit for some 
particular values oiuJi/uj2 & X CL K+ |FV| . but it does not assume vahdity of H17a|l 
in this regime. Our result consists thus of the observation that, after an appropri- 
ate rewriting, this modular transformation provides a representation for the elliptic 
gamma function that is well defined for all U)\li02 > |S2| . 



Below we shall need functional equations satisfied by S{u; uj) 

S{u + uji;uj) 1 S{u + LU2',<^) 1 



27ri- 



(19) 



e 



and its asymptotics for u going to infinity. Let us fix the incommensurate quasiperi- 
ods LOi and UJ2 such that Im(ti;i/ti;2) > (so |g| < 1). Then it follows from the infinite 
product representation (|16|l and its modular inverse 

( -2Triu/uJi . ~\ 

Siu;^,,U2) = e-"^-^-("^-) ;"_ (20) 

(e ^""/"2g;g)oo 

that 

lim S{u;uj)^l, (21a) 

lim e'^*-^^.^("''*')S'(M;a;) = 1. (21b) 
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It turns out that the same asymptotics also holds for the boundary domain uji /a>2 > 
(so |g| = 1), as can be verified by means of an integral representation for S{u] uj) 

IKtnilfea] . 

Remark 2. If we take the limit p,r —> in 117a|l . then the transition from G(u;uj) 
to the double sine function simplifies to 

T / —7TiLJ3 — — ^T^/ — 27rz— — 27r't— — 27ri— > 

lim le 121^11^2 1 (e "3 ; e , e "2 

Im(^),Im(^) — + 00 V 

^ ^-^^'''"'^'^Zllr'^^' S-\u;lui,uj2). (22) 

Such a limiting relation was first derived in a different way and in a stronger sense 
by Ruijsenaars |Rul |. 

4. The elliptic beta integral on the unit circle 



We now turn to the elliptic beta integrals. Let \q\, \p\ < 1 and let i„, n = 0, . . . , 4, 
be five complex parameters satisfying the inequalities |t„| < 1 and \pq\ < \A\, where 
A = Y\n=Q ^n- The elliptic beta integral of [HI] states that 

1 /■ n^=or(in^'^;g,p) _ '^U.0<7i<m<4^itntm;q,p) ^^^^ 



2m Jr r(z±2,Az±;g,p) z (g; <z)oo(p;p)oo DLo r(^in^ 9,p) ' 

where T denotes the positively oriented unit circle. Here we have employed the 
shorthand notations 

m 

r(zi, . . .,z^;q,p) =Y[^{zi;q,p), 

1=1 

r{tz'^;q,p) = r{tz,tz~'^;q,p), r(z=^^;q,p) ee r{z^ , z~^; q,p). 

For p = 0, the integration formula H23|l amounts to an integral explicitly constructed 
by Rahman in |R,ah| through a specialization of the Nassrallah-Rahman integral 
representation for a very-well-poised basic hypergeometric 8<^7 series. 

The following theorem provides a modified elliptic beta integral that — unlike 
— is well defined for |g| = 1. 

Theorem 2. Let Im(uji/LU2) > and Im[tjj^ / loi) , Im{Lj'i / UJ2) > 0, and let gn, n = 
0, . . . ,4, be five complex parameters subject to the constraints 

Im{gn/uj3) < 0, Im{{A - uJi - ^^2)7^3) > 0, 

with A = X]n=o 3"- Then the following integration formula holds 

n<m<4 

G{g 



-u,i2 G{±2u,A±u-^) LJ2 nLoG(^-5n;^) 
where 

_ -2(g;g)oo 



(24a) 



(24b) 



(9;9)oo(p;p)oo(r-;r)oo 

Here the integration is meant along the straight line segment connecting —u-^jl to 
L0-il2 and we employed the shorthand notation G{a ± 6; a;) = G{a + b,a — b]Oj). 
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Proof. We start by substituting relation ()17a|l into the left-hand side of H24a|l . This 
yields 



where 



\ n=0 / \ n=0 



(25) 



Vm— 1 



The constraints on the parameters permit to employ formula l|23|) with the substi- 
tutions 

2; ^ e , i„ — > e , p e "3 ^ q ^ e "3 , 

which yields for (|25|l 



2^3^2~^e"° no<n<m<4r(e '^"'"-3'';f,p) 
(r~; f )^ (p; p)^ j^4^^ r(e-'"*^ ; f , p) 

no<ri<m<4 

G{g 

nLoG(^-5«;'^) ^ 



with 



t^2(r;f)oo(p;p)oo 

After applying modular transformation Q to the infinite products appearing in k, 
we obtain 

V «W2 (r; r)oo(p;p) 00 

One more application of (jJI) allows us to replace the exponential function by a ratio 
of infinite products, which entails the desired form of n in (j24bp . □ 



Let us now consider the formal limit p, r ^ of the integral (|24a|l . To this end, 
we fix the quasiperiods wi^2 such that Im(wi/cL'2) > and Re(cL'i/w2) > 0, and we 
furthermore pick = ituj2 with i > 0. For t +00 the integral of Theorem [21 
then degenerates formally to 

S{±2u,A±u;ui) du ^ o (g;g)oo Ut^oSjA- gn^t^) ^^5) 

n<m<4 

where the integration is along the line L = iuj2^, and with parameters subject 
to the constraints Re(g„/u;2) > and Re((.4 — L0i)/ui2) < 1- This integral was 
deduced by a similar formal limit from the standard elliptic beta integral (|23|l and 
rigorously proved by Stokman in (S^, where it was referred to as the hyperbolic 
Nassrallah-Rahman integral. 

Remark 3. The generalized gamma function [z;t]oo used in St coincides with the 
double sine function IjlGI) upon setting z ~ {lo2 — u)/uj2 and t = —102/^1- 
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Remark 4. To further elucidate the intimate relation between the elliptic beta 
integrals (|28|l and l|24a|l . let us recall that — according to the general definition 
introduced in llS2| — a contour integral is called an elliptic hypergeometric integral 
if, for some displacement loi, the ratio of its integrands A(w + aJi)/A(u) constitutes 
an elliptic function of u with periods 102, (say). Now, by the change of variables 

277-1 2-Ki 2ni a 

z — e^2 ^ t,-^ — e , A — ^ 
the integral H23() passes into the additive form 



^W— = = „ -■4-3„ : (27a) 

-.2/2 ^2 (g;q)„o(p;p)oon„=or(e'"~;<?,p) 



with the integrand given by 

T(e "2 ,e -2 ;q^p) 
where r(e°^''; 9,p) = ^(e"+^ e°"^ q,p). We then have that 

A(u+_wi) ^ 27r»i^^ 6'(e "2 ;p)g(e "2 ;p) -TT 6>(e -2 ;p) 

6'(e -2;p)6'(e -2 ;p) „=o 6'(e -2 ;p) 

which is seen to be an elliptic function of u with the periods 0^2 and uj-j,. With the aid 
of the difference equation 1)1 4a|) for the modified gamma function, it is not difficult 
to verify that the integrand of modified elliptic beta integral (|24a|l also provides a 
solution to (|28|l . Hence both integrals (|23|l and H24al) are elliptic hypergeometric 
integrals with integrands satisfying H28|) . Whereas the solution of 128|l originating 
from integral 123|l is well-defined only for jql < 1 (or for j^j > 1 upon performing 
the inversion q ^ q^^ \vl the elliptic gamma function jS2j ). the modified integral 
(|24ajl corresponds to a solution that extends analytically from the regime |g| < 1 
to the unit circle j^l = 1. In this sense the modified integral H24a|l may be seen as 
a |g| = 1 analog of the original elliptic beta integral 

5. Multiple integrals 

The following integral constitutes a multidimensional generalization of the ellip- 
tic beta integral in 



1 /■ T-r r(tef z^;g,p) ^ Jln^o r(tnzf ; g,p) dz^ dz_ 



rf-y^-y^- „ 11 



AT 



(2^z)^ Vi<^^t<^ r(z±z±;g,p) r(zf ,i?z±;g,p) zi z^ 

2^m T{P;q,p) Uo<n<.r<4nt'-'tntrn;q,p) 



where B = t^^~^ Ylt.=o ^(te^a;^; g,p) = T{tzx,tzx~^ ,tz^^x, tz~^x~^;q,p), and 
with parameters subject to the constraints \p\, \q\, \t\, |t„| < 1 and \pq\ < \B\. This 
multiple elliptic beta integral was first formulated as a conjecture in |DS1| . Next, 
it was shown in |DS2| that the conjecture in question follows from a vanishing 
hypothesis for a related multiparameter elliptic beta integral. Recently, a complete 
proof of the integral in was found by Rains 0. 
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The following theorem provides the corresponding multidimensional generaliza- 
tion of the modified elliptic beta integral in Theorem El 

Theorem 3. Let Im[L0i/uj2) > and Im{LU3/uJi), Im{uJ3/LU2) > 0, and let g, gn, 
rt = 0, . . . ,4, be six complex parameters subject to the constraints 

Im{g/uj3), Im{g,Juj3) < 0, Im{{B - Wi - ^2)7^^3) > 0, 

with B = {2N - 2)g + ^^^^ g,, . Then 

TT G(.g ± Uj ± Uk] (Jj) ^ n«=0 G{gn ± Uj]Uj) dui dujq 



Proof. The proof is analogous to that of Theorem |21 Specifically, after substituting 
(|17af) into the left-hand side of H30() and application of the multiple beta integral in 
(|29|l . we arrive at the integration formula stated in the theorem upon expressing the 
resulting evaluation in terms of the modified gamma function G{u] uj). To infer the 
correctness of the value of the proportionality constant k^N\, one observes that the 
dependence on g in the factors originating from the exponential multipliers cancels 
out. The expression for the proportionality constant then follows from the fact 
that for 5-^0 integral ijSO)! reduces to the A''-th power of the elliptic beta integral 



Remark 5 . In |DS2[ IS2| various other types of multiple elliptic beta integrals were 
formulated. These can be extended to the unit circle \q\ = 1 in a similar fashion. 

For p = 0, elliptic beta integral reduces to a Gustafson's multiple integral 
corresponding to the Nassrallah-Rahman type generalization of the Selberg integral 
jGu2| . The following theorem provides a corresponding multiple analog of the 
integral in 1)26(1 . The integration formula in question can be obtained formally from 
the modified elliptic beta integral H3U|) with LOxjioi > by taking the limit p,r ^ 
in the manner explained below Theorem [2| 

Theorem 4. Let 101,102 be quasiperiods such that Im{uji / L02) > and Re{LOi / UJ2) > 
0. Furthermore, let g,gn, n = 0, . . . ,4, be parameters subject to the constraints 
Re{g/LOi), Re{g/L02), Re{gn/L02) > and Re({B~bJi) / UJ2) < 1 (with B as in Theorem 
0). Then 





(1^ 



□ 




(31a) 



where L — 1102^^ 



A(u;g) 



n 



^(ibuj ± Ufc; w) Yj S{±2uj, B ± Uj'jU:) 



(31b) 



l<j<k<N 



S{g ± u, ±uu;i^) nLo S{gn ± u,;u) 



and 



(g; g) 

(9; 9) 



N N 



SUa;^) no<n<m<4'S'((j - 1)5 + 5™ + ffm;w)' 



A/'(g) = (-2) 




(31c) 
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Through a parameter speciahzation, Gustafson's muhiplc integral of the Nassral- 
lah-Rahman type can be reduced to a multiple Askey- Wilson integral first derived 
in |Gulj . The corresponding degeneration of Theorem 0] reads as follows. 

Theorem 5. Let cui, lu2 be quasiperiods such that Im{LOi/ LJ2) > and Re{LOi/ LJ2) > 
0, and let g,gn, n = 0, . . . , 3, be parameters subject to the constraints Re{g/uji), 

Reig/i02), Re{gn/i02) >0 and Re{{B - UJ2) / u^i) < 1 with B = {2N -2)g + Y.l=ogn- 
Then 

T-r 5(±Uj ± life; u)) 1^ S{±2uj]u}) dui du^ 

i<}-}k<N s{g±u,±u,;u,) nLo5(.9n±^.;'^) 

= ( o'i^Ari (g'g)^ TT -^(ff'^) S{{1- i)g + B;uj) 

For = 1, the integral in Theorem [S] reduces to a single variable Askey- Wilson 
type integral 

S{±2u]u>) du _ ^ {q;q)oo S {gp + gi + g2 + gs; ^gg^ 



n<m<3 

Sig n I gm 7 ^ ) 



which was established by Ruijsenaars jRu3| and Stokman jSt| . 

Formally, the integral of Theorem|31follows from that of Theorem0]with uJi/u!2 > 
upon setting — + (74 + 11^2^ and performing the limit t — + +00. However, it is not 
a simple matter to upgrade such formal limiting relations between Theorem |21 and 
Theorems0I|Slto rigorous proofs of the latter integration formulas. A direct proof of 
Theorems 0] and O modelled after Rains' proof of the multiple elliptic beta integral 
l)29|l . is given in Section]^ below. As it was communicated to us by Stokman after 
finishing this paper, Theorem can be proved by a multivariable generalization of 
the method of Stj as well. 



6. Proof of Theorems 01 and [51 

We first detail the proof of Theorem 01 and then indicate some modifications so 
as to incorporate Theorem O 

Let us for the moment assume that the quasiperiods Wi, 0^2 are incommensurate 
over Q. The double sine function S{u; lj) then has simple zeros located at the 
points u = —ujiN — W2N and simple poles at m = ti)i(l + N) + 0^2(1 + N). Therefore, 
the integrand A(u; g) in H31a|) has poles at the points 

±uj = -B + LUi{l + N)+LU2{l + N), gn + LOiN + UJ2N, n = 0,...,4, 
g + Uk + wiN + g - life + cJiN + W2N, A: = 1, . . . , A^, fc ^ j, 

where j = 1, . . . , A^. 

Combination with the asymptotics in Eqs. (|21a|l . (|21b|l reveals that the quotients 
S{u\ u})/S{g + u; uj) and S{2u, B + u; i^)/Y[n=o S{9n + u) are smooth and bounded 
on the complex line L = 1^2]^. Indeed, for u — ixuj2, a; e M we stay away from 
poles and we have that 

5(1^22;;'*^) J 0(1) for .T ^ +00 

S{g + iuj2x;u:) ~ \ 0(5-2 W'^i ^ for a; ^ -00 
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and 

S {2iijj2X , B + iijJ2x; ijj) I 0(1) for a; ^ +00 



nLo^(3"+*^2x;c^) 1o(e2-(2(^-i)9M+i+-.M)) for 



X —00 



It thus follows that the integrand A(u;g) is smooth and exponentially decaying at 
infinity on the integration domain L^. Hence, the integral in Eq. H31a|) converges. 

To infer the validity of the integration formula we distinguish three parameters 
9o,gi,92 and factorize the integrand as A(u;g) = A+(u)A_(u) with 

A fu) - TT -^("j ^ "fc;^) TT S{2uj,B - Uj,uji + C-Uj;u}) 

" i<j<t<N + u, ± ufc; u:) n s{u^ +C + uf,u) HLo S[gn + u.-u,) ' 

(34) 

where C — {N — l).g + .go +51+32 and A_ (u) = A+(— ui, . . . , — it„). Similarly, we 
introduce the shifted functions 

n (35) 

JX- s{2uJ,B + ^^-UJ,c + ^^~uf,u}) 

^ 11 



and A_(u) — A_|_(— ui, . . . , — u„), which provide a factorization of the integrand 
for a shifted set of parameter values: 

A / \ A / \ A / ^1 ^1 ^1 ^1 '^l ^ 

^+(u)^-(u) = A(u;.g,5ro + Y,ffi + Y,ff2 + y^93 - Y'^i - y)- 

Let us for the moment assume that the parameters are such that the shifted pa- 
rameters (as well as the ones obtained after shifts by ±012/2) also belong to the 
domain indicated in the theorem. We then have the following equality 

/ A+{ui + '^, . . . ,UN + ^)A_(u)dui • ■ -duN 

= / A+{u)A^{ui - . . . ,UN - '^)dui - ■ - duN, (36) 

which follows by shifting the integration contours L on the left-hand side by —uji/2. 
Notice that such shifts are permitted by the Cauchy theorem due to the absence 
of poles in the strip between L and L — a>i/2 combined with the exponential decay 
at infinity. Indeed, the quotient S{u;uj)/S{g + u;u)) is holomorphic on the strip 
{u — sui + ixu!2 I < s < 1, —00 < X < 00} and the quotient S{2u,B + 

u,C + u;lj)/ (5(C + ^+u,g3-^ + u,g4-^+u;uj) JlLo ^(.9n + ^ + ^)) 
is holomorphic on the strip {u = slui + ixuj2 |0<s<l/2, — (X)<a;<oo}. 
By performing sign flips of the form Uj —Uj in the integration variables and 
summing over all 2^ possible ways, we obtain from (|36|1 that 

/ p(u;g)A+(u)A_(u)dMi---duAr 

= / p{u;g)A+{u)A^{u)dui---duN, (37) 
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with 

K+{v^Ul + ^^,...,v^uN + ^t) 



p(u;g)= XI 



£=l,.-.,Af 



p(u;g)= 2^ - 



^ = l,...,Af 



2 

A_(l/ilti,...,r/ArMjv) 



Simplification of p(u; g) reveals that 



- 1 _ 

■'f = ±i l<i<fe<7V ^fe j = l 

1=1, -...N ~ ~ 

N 

= \{{i-t^-hQh)(i-t^-h^,t2)(i-t^-hit2), (38) 

with i e^'^^f/"^ t„ = e^'^^S"/"^^ Zk = e^'^™''/'^^ The summand in ^ is invariant 
under permutations of Zk and inversions Zk ^ \ ^ ■ The product of p(u; g) and the 
factor 

-Q (1 - ZjZk){l - ZjZ^^) -Q ^ " 

l<j<k<N 3 = 1 '^^ 

yields a Laurent polynomial in Zj, j — 1, . . . , TV, which is antisymmetric with respect 
to both transformations (separate permutations of zj and inversions zj — > zj^). 
Any such polynomial is proportional to the multiplicative factor given above. The 
constant of proportionality is found after setting Zj — tot^^^, which leaves only 
one term in the sum (with all vj — 1) giving the right-hand side expression in (|38|l . 
In the same way, one obtains for p(u;g) that replacing in relation (|38() to, 1,2 by 

J- „-l/2 iN-l-i J- ^ „l/2 
13,4(7 ' , t tOtlE2'7 ' J 

. 5: n ^n( "-'-'7"i; 

-f = ±i l<j<fe<7V i fc 3 = 1 \ ^ ^3 

X (1 - t^-H,t^Wl^z^^){l Sg-i/^z-^O^ 

N 

where B = e^^^^l"^^ . The derived expressions demonstrate that the functions in 
question are constant in the integration variables u, that is p(u; g) = p(g) and 
p(u; g) — /5(g). Hence, we can pull the corresponding factors out of the integrals 
and rewrite H37|) as 

/ A+(u)A_(u)dwi---duAr = 4^ / A+(u)A_(u)dMi---duAr, 
A" P(g) A" 

whence 

AA(g^ ^ p(g) 

A^(5,5o + ^,<?i + f,52 + f ,33-^,54-^) p{&)' ^ ' 
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As a result, we deduce that the ratio of the left- and right-hand sides of (|31a|l is 
invariant with respect to the shifts 50,1,2 <?o, 1,2+^1/2, (73,4 53,4 — ^1/2, and, by 
symmetry, any permutation of indices of the parameters. The double sine function 
S{u; uj) and, so, the integrand in (|31a|l and the integral's value A/'(g) are symmetric 
with respect to the permutation of uji and uj2 |KLSj . The contour of integration 
L breaks the symmetry between uji.2, but the transformations used in H34|l - 139() 
are purely algebraic and do not depend on the contour of integration. Therefore, 
the ratio of interest is invariant with respect to the shifts 50,1.2 ~* 3q.i,2 + '-^2/2, 
53,4 ~* 53,4 ~ ^2/2 as well (with all permutations of indices of parameters). 

By analyticity, without changing the integral's value we can replace the contour 
of integration L by any other contour embracing the same set of poles. For an 
appropriately deformed contour, we can establish invariance of the ratio of interest 
under the shifts gn — > gn + kLUi/2 + muj2/'2, for arbitrary k,m gZ. Taking uji,lu 2 > 
0, we can choose a subset of these points with a limiting point in the parameter space 
for which we can choose L as the integration contour. Moreover, making sequential 
^1^2/2 shifts in different directions we can escape large intermediate deformations 
of the integration contour (such an argument is similar to the one given in ^SIJ. 
Therefore, the ratio of the left- and right-hand sides of (|31a|l is equal to a function 
of and g, which we denote as f{uji,uj2,g)- 

In order to see that f{uji,ijj2,g) actually equals to one, it is necessary to use 
an analog of the residue formula derived in |US1| . Namely, we take one of the 
parameters, say, go such that one pole from the half plane Re(u/w2) > crosses 
the contour of integration L. A similar move takes place for the pole at m = —go 
from the Re(u/w2) < half plane (due to the reflection invariance). We keep intact 
all other poles in the half planes to the left or right of L. This is possible to do 
by an appropriate choice of 5, gi, . . . , (74, and loi.2- It is not difhcult to see that the 
residues of these crossing poles taken, say, over the variable u n have poles at the 
points Uk = ±(50 + g): k = 1, ■ ■ ■ , N — 1, (instead of Uk = ±go) and they are still 
located to the left or right of L for Re(((7o + g) /UJ2) < 0. Similar shifts of the poles 
occur each time we calculate the residues. Therefore we denote pk = go + (fc — l)g 
and impose the restrictions 



Because of the taken constraints upon g, 101^2, we get a simpler residue formula 
than the one derived in |DSlj : 



where is a deformation of the contour L such that it separates the same sets of 
poles as L did before we started to change go- The factor 2™ emerges because the 
residues appear in pairs and their values coincide (due to the Uk —>■ —Uk reflection 
invariance of the integrand and different orientation of the contours encircling poles 



of m cycles and the number of ways to pick up these cycles out of N possibilities. 





(40) 




count the number of orderings 
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The residue functions have the form /io(u) = A(u; g) and for m > 

P„i{u) = K„(5,„^Ar_„(u)A7V-m(u;g), (41) 

where Ajv_„j(u; g) is obtained from integrand Ij31b|l if we replace in it by iV — to 
but keep B = [2N — 2)g + X]n=o 3" unchanged. Other coefRcients are 

^( ,Nm (g;g)S TT S{±pk-Pj\i^) ^ S{~2pi,B±pi\U}) 

('^;'^)S,</<L^(5±P.-P.;-)finLi^(5«±P.;-)' 

and 

A r,,^ - TT S{±pj±uk;i^) 

Om.N-mW - I I „, I — T- 

ih<m S{g±p,±u,;u:) 

l<k<N-m 

The expressions for /im(u) are derived by induction. Indeed, the form of pi{u.) is 
easily established after taking into account the relation 

lim " T .90 _ _^ (<?; g)oo 



and the fact that the contours encircling the corresponding poles are oriented clock- 
wise for the upper signs and anticlockwise for the lower signs (this gives the total 
minus sign in 

Suppose that p,n is given by (|41|) for some to > 1. In order to find Pm+i it is 
necessary to compute the residues for poles located at ujv-m = ipm+i- A simple 
computation shows that, indeed. 



where Cm is a small size clockwise orientated contour encircling the pole at uty— rn — 

Pm+l- 

By analyticity, our deformations of the parameters and of the contour of inte- 
gration do not change the integral value and, therefore, the right-hand side sum in 
(En|l equals to /(wi, cj2, g)A/'(g). We now divide both sides of this equality by 7V(g) 
and take the limit 54 ^ —g^ — {N — l)g. For m < N, the coefficients /tm(g), which 
can be represented in the form 

_ (_-. ym (g; g)S TT / S {g , {2 - m - l)g - 2go; uj) 

S{2go + Y.t=i gr + {2n + l- 3)g, E^^i 9r + {2n - I - l),g; uj) \ 

^ 4 I ' 

Ylr=iS{9r +go + {l- 1)9, 9r ~ 9o ~ {I ~ ^)9i^) J 

do not contain diverging factors in this limit and the integrals, which they are 
multiplied by, remain bounded. Therefore, only the term with m = N survives 
and, by simple computation, we obtain 

lim — , = 1, 

34— 50-(W-l)g A/(g) 

which means that /(wi, 0^2,5) — 1- After proving equality H31a|) in the taken 
restricted region of parameters (where the parameters shifted by ±a;i_2/2 satisfy 
the needed constraints and wi_2 > 0), we can analytically extend it to the values 
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of U!i^2 and parameters g,gn in the domain indicated in the formulation of the 
theorem. Theorem 0] is thus proved. 

We now turn to the Askey- Wilson type integral (|31a|l . Its convergence conditions 
essentially differ from the previous case. Indeed, we have 

S{2iijj2X] (jj) J 0(1) for X +oo 

Y{l=oS{9n + ioJ2X\u:) ~ jo(e2--(i+"^/'^i-S'=off"/-i)) for a; ^ -oo ' 

Combining together these limiting relations with the asymptotics for the ratio 
S{iuj2X]oj) / S{g + iuj2x;u3), we see that the integrand remains bounded in the 
integration domain and decays exponentially fast on its infinities if we take 
Re((B - L02)luji) < 1, where B ^ {2N - 2)g + ELo5«- 

Invariance of the ratio of left- and right-hand sides of equality H31a() under the 
specified parameter shifts relied only on algebraic manipulations with the integrand. 
Therefore we can repeat them for the limiting expression of the integrand appearing 
after taking the limit Im((74/a-'2), Im((74/wi) — > -|-oo (or — > 0). This simplifies 
the integrand for (|31a|) to the one for (|32|l . Therefore, limiting analogs of equalities 
show that the ratio of the left- and right-hand sides of do not depend 
on the shifts in the parameter space 50.1,2 ~* 3o,i,2+<^i,2/2, 53 — > 53—^1^2/2 and the 
ones obtained by permutation of indices. Using, again, an analytical continuation 
and the appropriately simplified version of the residue calculus, we see that the ratio 
of interest is actually equal to one. As a result, we establish validity of integral 1321) 
as well. 
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